Coarse grained dynamics of the freely cooling granular gas in one dimension 
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We study the dynamics and structure of clusters in the inhomogeneous clustered regime of a freely 
cooling granular gas of point particles in one dimension. The coefficient of restitution is modeled as 
ro < 1 or 1 depending on whether the relative speed is greater or smaller than a velocity scale 5. The 
effective fragmentation rate of a cluster is shown to rise sharply beyond a 5 dependent time scale. 
This crossover is coincident with the velocity fluctuations within a cluster becoming order 5. Beyond 
this crossover time, the cluster size distribution develops a nontrivial power law distribution, whose 
scaling properties are related to those of the velocity fluctuations. We argue that these underlying 
features are responsible behind the recently observed nontrivial coarsening behaviour in the one 
dimensional freely cooling granular gas. 

PACS numbers: 45.70.Mg, 45.70.Qj, 05.70.Ln 



I. INTRODUCTION 



Consider a collection of particles, initially distributed 
randomly in space, evolving in time through ballistic 
transport and inelastic collisions. Such a system has been 
studied extensively as a simple model of granular systems 
as well as a tractable model in nonequilibrium statistical 
mechanics [l|-[T2|. At initial times, the system undergoes 
homogeneous cooling, the energy decreasing with time 
t as t~'^ in accordance with HafFs law [l^ |l3]- In this 
regime, the particles remain homogeneously distributed 
with inter-particle spacing being the only relevant length 
scale. At times larger than a crossover time tc, the sys- 
tem crosses over to an inhomogeneous clustering regime 
where the energy decreases as , where 9 varies with di- 
mension and is different from 2 in dimensions lower than 
the upper critical dimension. In this regime, there is a 
growing length scale £<, determined by the size of the 
largest cluster. 

In one dimension, much more is known than in higher 
dimensions. Through an exact solution of the prob- 
lem with coefficient of restitution set to zero (sticky gas), 
and extensive simulations 6] of the inelastic gas, it is 
known that 9 — 2/3. The sticky limit may also be 
mapped to the dynamics of shocks in the inviscid Burg- 
ers equation [l^. Recently, we showed that, when the 
coefficient of restitution depends on the impact velocity, 
then a new timescale ti further subdivides the inhomo- 
geneous clustering regime into two sub- regimes [ill, ■ 
This was based on a study of the density-density and 



•Electronic address; |mahendra.statp @gmail.com| 
t Electronic address: dibyenduQphy.ii tb.ac.in'] 
•f Electronic address: ,rrajesh@imsc.res.in| 



the velocity- velocity correlation functions. For times 
tc < t < ti, these structure functions scale exactly as 
in the sticky gas, obeying what is known as the Porod 
law [l^. However, for times t > ti, the inelastic gas 
deviates from the sticky gas limit and the correlation 
functions violate Porod law. In addition, the density dis- 
tribution and inter-particle distance distribution develop 
into power laws that are qualitatively different from that 
seen in earlier times. We will refer to the two sub-regimes 
as the Porod and fluctuation dominated ordering regimes, 
respectively. 

Although the macroscopic statistical quantities studied 
in Refs. [HI, [13] establish that the Porod and the fluctua- 
tion dominated ordering regimes are distinct, they do not 
reveal how fluctuations start dominating beyond the time 
scale ti . It was speculated that the fluctuation dominated 
ordering regime should have an effective process of frag- 
mentation that will compete with the otherwise strong 
effective aggregation. Whether coarse grained density 
clusters break up or remain coherently moving objects 
can be directly checked by studying their dynamics. In 
this paper, we study cluster dynamics, in particular ef- 
fective fragmentation rates, and show that the ordering 
process gets disturbed beyond a certain timescale. 

The second question is regarding the origin of the 
crossover timescale ti . Earlier this was shown to depend 
on a velocity scale S associated with the coefficient of 
restitution [ill, [i3| • The coefficient of restitution is often 
modelled as a function of the relative velocity- rather 
than being a constant- such that collisions become near 
elastic for relative velocities smaller than S. This prop- 
erty is consistent with experiments [T7l - [l9j . required by 
theory as well as essential in simulations to pre- 

vent inelastic collapse [2|. In this paper we demonstrate 
that the time ti is marked by the particle velocity fluc- 
tuations within density clusters becoming of the order 6. 



2 



Thus, it is not the typical velocities but rather the veloc- 
ity fluctuations that matter for correlation functions. We 
also present a consistent scaling theory to understand the 
dependence of velocity fluctuations, fragmentation rates 
and cluster size distribution on the parameters S, time t 
and cluster size m. 

In Sec. ini we define the model, the quantities of in- 
terest, and give details of the simulation. In Sec. IIIIl we 
present results from numerical simulations for velocity 
fluctuations, fragmentation rates and cluster size distri- 
bution. We develop a scaling theory which enables us 
to understanding the scaling of the above quantities in 
terms of two exponents. Sec. II VI contains a summary and 
discussion of results. 



II. MODEL AND DEFINITIONS 

In this section, we define the microscopic model, and 
the quantities of interest. The model consists of a collec- 
tion of N point particles of equal mass on a ring of length 
L. Each particle moves ballistically until it collides with 
another particle. The collisions conserve momentum, but 
are inelastic, such that when two particles with initial ve- 
locities Ui and Uj collide, the final velocities u[ and u'j are 



determined by 



1-r 



1 + r 



(1) 



where r is the coefficient of restitution. The coefficient of 
restitution is velocity dependent such that collisions be- 
come elastic when the relative velocity tends to zero. As 
mentioned in the introduction, this feature is consistent 
with experiments and theory, and often used in simula- 
tions to circumvent inelastic collapse. Our results will 
not depend on the detailed dependence of the coefficient 
of restitution on relative velocity: we therefore choose a 
convenient functional form 0, 01, 



ro if Vrd > S, 

1 if Vrc\ < S, 



(2) 



where (5 is a velocity scale in the problem. The colli- 
sions are elastic for relative velocities smaller than 6. The 
particles are initially distributed randomly in space with 
their velocities drawn from a Gaussian distribution. 

Starting from the above microscopic model, we desire 
to study emergent processes such as fragmentation and 
aggregation of a collection of particles. To this end, a 
coarse grained description has to be introduced, which we 
do as follows. Divide the ring into N equally sized boxes. 



Let the number of particles in the i 



th 



1,2,.. 



,N) 



box be called the box density. We define a cluster to 
be a collection of contiguous boxes with non-zero box 
density surrounded by two empty boxes. The total num- 
ber of particles in this cluster will be called the size of 
the cluster. Similar definition has been used elsewhere, 



for example see !21l|. In earlier papers |l2|, we had 
studied velocity-velocity and density-density correlations 
using the coarse grained box density. However, a single 
particle moving across the boundary of a box results in 
change of density. To study clusters, any coarse graining 
should make sure that particles move significantly, i.e. 
by at least a box spacing [L/N) to cause a change in 
configuration. The above definition of a cluster has this 
property. 

One of the quantities of interest in this paper is veloc- 
ity fluctuations a within a cluster. Let Uc be the centre 
of mass velocity of a cluster, i.e., Uc = 'm~^Y^^l'^i^ 
where are the velocities of the particles constituting 
the cluster, and m is the size of the cluster. Then, 
the velocity fluctuations of that cluster is cr(m, t) — 
~ Uc)^, where the summation is again over 
all the particles constituting the cluster. When a is much 
smaller than m^, then the cluster is compact and stable 
with respect to the velocity fluctuations. However, if they 
are of same order, or ii a ^ 8"^ , then the cluster may start 
breaking apart. This effect can be captured by defining 
an effective fragmentation rate of a cluster, as described 
below. 

In the space of cluster sizes, a stochastic dynamics may 
be defined, with rates for the different processes being 
determined from simulations. As the clusters evolve in 
time due to the entry and exit of particles into and out 
of a cluster, we ask what the rates of transition from a 
cluster size m to m! are. Let W{m — >■ m']t) be the rate 
at which a size m changes into a size m' at time t. If 
N{m, t) is the number of clusters per unit lattice site of 
size m at time t, then its time evolution is described by 
an effective Master equation ^22]: 

'^^^^'^^ - ^M^(TO'^m;t)iV(TO',t) 

- ^Wim-^ m';t)N{m,t). (3) 



If in a time interval At, the cluster size decreases, 
then the cluster is said to have undergone fragmenta- 
tion. Thus, an effective fragmentation rate Wf{m,t) of 
a cluster of size size m at time t can be defined as 



Wf{m,t) = W{m 



;t). 



(4) 



m' <m 



We note that, in the above, fragmentation is an emergent 
process, not defined apriori in the microscopic dynamics, 
unlike some other models of gran ular gas where fragmen- 



tation occurs on collision [23|, |24| 



We study the model by means of event driven molec- 
ular dynamics simulations p5l |. In the simulations the 
number density N/ L is fixed to be to one and the num- 
ber of particles to be iV = 20000. The results in the paper 
do not depend on the precise values of the parameters rg 
and 6, as long as 6 is much smaller than initial velocity 
differences of adjacent particles. We use generic values 
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6 = 0.001, 0.002, 0.004, 0.008, and ro = 0.1 in the simula- 
tions. The initial velocities are chosen from a Gaussian 
distribution with width 1. The data is typically averaged 
over 20000 - 30000 different initial conditions. All aver- 
ages will be over space and different histories and will be 
denoted by (• • • ) . Also, we use reduced units in which all 
lengths are measured in terms of initial mean inter par- 
ticle spacing and times in terms of initial mean collision 
time. 



III. RESULTS 

There are four velocity scales in the problem. First is 
the typical speed of a cluster which decreases in time as 
Second is the root mean square velocity fluc- 
tuations y/a within a cluster (discussed in detail below) . 
Third is 6, characterising the coefficient of restitution [see 
Eq. ([2])], while the fourth corresponds to the initial ve- 
locity distribution. At large times, there is no memory of 
the initial velocity distribution, and it will play no role 
in the subsequent discussion. When the typical speeds 
become of order 5, then almost all collisions are elastic 
and energy does not decrease any more. We will denote 
the latter crossover time by t2- Clearly, t2 ^ S~^. It 
is possible that the velocity fluctuations scale with time 
differently from the typical velocity. If so, we have a pos- 
sibility of a different crossover time which is marked by 
the velocity fluctuations becoming order 6. 

We first characterise the velocity fluctuations a. In 
Fig. [1] we show the variation of a with time t for dif- 
ferent values of S and two values of cluster size m. For 
short times, a is independent of S, and after initial tran- 
sients, decays in time as a power law, with the exponent 
independent of m and the prefactor dependent on m. At 
large times, a deviates from the power law behaviour and 
is constant for a while. We argue that this crossover oc- 
curs when a is of order — velocity fluctuations and 
hence relative velocities are such that collisions within 
a cluster become near elastic. Elastic collisions tend to 
smoothen out density inhomogeneities. Therefore, clus- 
ters become less compact and fragmentation is initiated. 
The observations are mathematically summarised as 

a{t,S) ~ (5^1 (i'^^^'^') , fixed m, (5) 
a{t,m) ~ /2 , fixed (5, (6) 

where xi, X2 are scaling exponents and /i, /2 scaling 
functions such that fi{z) ^ z^^^, z <C 1, and f2{z) ^ 
z <C 1. Thus, for fixed 6, cr ^ m^^t~^^ for initial 

times. 

The exponents Xi, X2 may be obtained from the data 
collapse of the data in Fig. [T] when scaled as in Eqs. (O 
and dl]). The scaled data is shown in Fig. [^Ja) [Eq. ([5])] 
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FIG. 1: The velocity fluctuations a within a cluster of size m 
as a function of time t for difi'erent S. The data are for sizes 
m — 19 and m = 76. For short times the curves, for a given 
size, decay as a power law and is independent of 5. Deviation 
from power law is seen earlier for larger 5. 
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FIG. 2: Data coUapse when scaled velocity fluctuations a5~^ 
is plotted against (a) scaled time variable t5°'^^' and (b) scaled 
time variable fm""'**. The data are for (a) m — 19 and (b) 
S = 0.004. The solid straight lines have slope = —3.0. 

and Fig. [2ljb) [Eq. ([6])]. From these, we obtain 

xi = 3.00 ±0.06, (7a) 
X2 = 2.66 ±0.08 (7b) 

Note that these values of xi and X2 imply that the 
crossover time ti, relevant for cr, scales as S~'^^^^ ~ 
^-0.66 This time scale is much smaller that t2 ~ 5~^, 
which is the crossover time associated with velocities of 
nearly all particles becoming of order 6, i.e. all collisions 
becoming near elastic. 

To contrast the above intra-cluster velocity fluctua- 
tions a with typical cluster velocities, we study the centre 
of mass velocity of a cluster. It is well known that av- 
erage energy per particle in the cooling gas decreases as 
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FIG. 3: The square of the centre of mass velocity (u^) of a 
cluster of size m as a function of time t for different S. The 
data are for sizes m = 19 and m = 76. For the times shown, 
there is no dependence on 5. The solid line has slope —2/3. 
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FIG. 4: The velocity fluctuations amax of the largest cluster, 
as a function of time t for different S. For short times the 
curves are independent of 5. Deviation from power law is seen 
earlier for larger 5. The solid straight line has slope —1.22. 



t~^^^ 0. Not surprisingly, we find that the (it^) for a 
cluster decays with the same law. In Fig. [31 we show 
the variation of (u^) with time for different 5 and two 
different cluster sizes, (u^) decreases as t~^/^ at large 
t. There is no signature of any crossover across any in- 
termediate time scale ti nor any dependence on S. The 
intra-cluster near elastic collisions affect velocity fluctua- 
tions but not the typical speeds, which are affected only 
by cluster-cluster collisions. 

We provide a further check for the exponent values in 
Eq. (O by quantifying the velocity fluctuations (Tmax of 
the largest cluster in the system. From Eq. ^ , we obtain 
that CTmax ~ -^max^'^S whcrc Mynax{t) IS thc sizc of the 
largest cluster at time t. Noting that Mmax(i) ^ ^^^^ 
0,101, we obtain (Tmax ~ t'^^2/3-xi _ where we 

substituted the values of X2 and xi from Eqs. ^ and 
(|6]). In Fig. |4l we show the variation of Cmax with time 
t for different values of d. The temporal regime which is 
independent of S is consistent with the exponent 1.22. 

We now show that the above crossover of a from a 
power law is linked very closely to the initiation of frag- 
mentation in clusters. The fragmentation rate Wf{m,t) 
deflned in Eq. ^ is numerically measured as follows. At 
time t, all the clusters of a particular size m are identi- 
fled. At time t + At, the fraction of the identified clus- 
ters whose size has reduced is calculated. That fraction 
is equal to Wf{m, t)At. In the simulations, we choose At 
to be one, so that sufficient statistics may be obtained. 

In Fig. [5l we show the fragmentation rate Wf for clus- 
ter size 19 for different values of S. There is a sharp in- 
crease in the fragmentation rate, with the increase setting 
in earlier for larger S. In the inset of Fig. [Sj we superim- 
pose the fragmentation rate on the plot of a with time 
for the same value of d and cluster size. Clearly, the in- 
crease in fragmentation rate coincides with the deviation 
from power law behaviour of the velocity fiuctuations. 
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FIG. 5: Fragmentation rate Wf as a function of time t for 
different 5. The data are for cluster size 19. Inset: Veloc- 
ity fluctuations a and Wf as a function of time. Increase in 
fragmentation rate coincides with saturation of velocity fluc- 
tuations. The data are for S = 0.004 and cluster size 19. 

This increased fragmentation rate results in fiuctuation 
dominated coarsening and ultimately to the breakdown 
of Porod law [ll| . 

The fragmentation rate Wf depends on the cluster size 
too. In Fig.[6l we show the variation of Wf with t for dif- 
ferent m and fixed S. The crossover time beyond which 
increased fragmentation is seen increases with m. Know- 
ing that t ^ m^^/^i (Eq. ([6])), we write 

Wf (t, m) ~ m"/™ (^j-) , fixed S, (8) 

where rj is an unknown exponent. We see that for i] = 0.5, 
we obtain excellent data collapse [see inset of Fig. |B] . 

An increased fragmentation rate will result in clusters 
up to a size m*{t, 6) breaking apart. We believe that this 
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FIG. 6: Fragmentation rate as a function of time t for 
different cluster sizes m. The data are for 5 = 0.008. Inset: 
Data collapse when Wj and time are scaled as in Eq. ([S]) with 
r\ = 0.5. 
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FIG. 7: The average cluster distribution {N{m,t)) as a func- 
tion of cluster size m for (a) fixed S — 0.004, different times 
and (b) different S, fixed time t = 16384. 

is the origin of the breakdown of Porod law, as fragmen- 
tation results in new structures at small scales. One way 
to capture this is to study the average cluster size distri- 
bution {N(m, t)), where N{m, t) is the number of clusters 
of size m at time t and the average is over space and his- 
tories. We would like to investigate whether, even in the 
presence of fragmentation, the cluster size distribution 
can be described by the sticky gas. We argue below that 
while some regimes of {N{m,t)) resemble the sticky gas, 
other regimes differ, but their scaling may be obtained 
from that of a{m, t). The mean cluster size distributions 
are shown for different times in Fig. ^a) and for different 
S in Fig.ElJb). 

We note that for fixed time, {N{m, t)) for large masses 
has no dependence on S [see Fig. El^b)]. Thus, we ex- 
pect that for masses greater than a mass cutoff m*(i, S), 
fragmentation is not relevant and {N{m,t)) should have 



FIG. 8: Data collapse for large cluster sizes m when the clus- 
ter distribution {N{m, t)) and m are scaled as in Eq. ([9|. The 
data are for S = 0.004. 

the same scaling behaviour as in the sticky gas. For the 
sticky gas, it is known ^] that 

(^K0)-^/3(^), (9) 

where the scaling function f^iz) ^ z^^/^, z <C 1 and 
/3(z) — > for z ^ 1. For masses m > m*{t,6) for which 
fragmentation is not important, we confirm numerically 
that the same scaling holds. In Fig.[8l we scale the data of 
Fig. ^a) as in Eq. ^ and we see excellent data collapse 
for large cluster sizes, confirming that fragmentation can 
be neglected for cluster sizes larger than m*{t,6). 

We note that in Fig. |8l there is no data collapse for 
small cluster sizes, when the data is scaled as in Eq. ([9]), 
and thus small cluster sizes cannot be described by the 
sticky gas scaling. We argue that its scaling can be ob- 
tained from that of the scaling of a. For a fixed S, varying 
t, {N{m,t)) should have the scaling form 

1 / 777 \ 

{N{m,t)) ^ -U [j^) ' rn<^m*, fixed <5, (10) 

where a is an exponent which we determine by examining 
the large z behaviour of the scaling function fiiz). For 
large z, f4{z) should be such that it crosses over to the 
small z behaviour of /3(z). Thus /4(z) ~ for z » 1. 

Comparing the time dependence of Eqs. © and pUj) in 
the latter limit, we obtain 

a=l + — wl.56. (11) 

2X2 

The data when scaled as in Eqs. (|10p and dTTI) with xi 
and X2 as in Eq. ([7]) is shown in Fig. IH^a). We obtain 
good data collapse for the small cluster sizes, showing 
that knowing the scaling behaviour of a helps us ob- 
tain the scaling behaviour of {N{m,t)). The small z 
behaviour of /4(z) can be determined numerically. We 
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FIG. 9: Data collapse for small cluster sizes m when the clus- 
ter distribution {N{m,t)) and m are scaled as in (a) Eq. pop 
for fixed S = 0.004 and (b) Eq. ^ for fixed t = f6384. The 
solid line has slope f.75. 



find that f4{z) ^ z ^ with t — 1.75 ± 0.08 [see solid line 
in Fig.lija)]. 

The scaUng of the small cluster sizes with 5 can also be 
obtained from the scaling of a. Knowing that m* scales 
as (5^/^^, we write 

(7V(m, 5)) c (^) , m « m\ fixed t, (12) 

where the exponent/3 can be determined as above by con- 
straining the large z behaviour of the scaling function 
/5(z) to be the same as the small z behaviour of fd.{z). 
This immediately implies that 75(2;) ~ z"^/^ for 2; ^ 1 
and 



1 

X2 



0.37 



(13) 



The data for cluster size distribution when scaled as in 
Eqs. (fT2|) and (fTS]) with xi and X2 as in Eq. ([7]) is shown 
in Fig. injb) . We obtain reasonable data collapse for the 
small cluster sizes. However, given the range and quality 
of data, it is possible to obtain data collapse for a range 

of X2. 



IV. DISCUSSION 

To summarise, we studied velocity fluctuations and size 
distribution of clusters in a freely cooling granular gas 
in one dimensional ring evolving via ballistic motion and 
inelastic collisions. The coefficient of restitution was tq < 
1 for relative velocity greater than 5 and 1 otherwise. The 
aim of the paper was to understand the consequences of 
a non-zero 5 on the structure of clusters for large times. 

For granular gases with realistic velocity dependent 
coefficient of restitution, it was recently shown [ll|, [T^ 
that the nature of coarsening in the inhomogeneous cool- 
ing regime is not the same at all times. Beyond some 



crossover scale ti, the coarsening behaviour changes at 
the macroscopic level from one that obeys Porod law to 
one that violates Porod law. These interesting numeri- 
cal findings lacked a mesoscopic explanation of how and 
why the crossover occurs. The current paper provides 
an explanation. We demonstrate in this paper, that the 
transition from sticky gas regime to fluctuation domi- 
nated ordering regime within the inhomogeneous cool- 
ing regime may be viewed as a growing dominance of an 
underlying fragmentation process competing against the 
dominant clustering process. The fact that clusters break 
up is shown by the change of behaviour of the variance 
of particle velocities within a cluster. This crossover in 
velocity fluctuations coincide with an increase in the frag- 
mentation rate of clusters leading to a richer fine struc- 
ture reflected in the density-density correlations. 

The velocity fluctuations within a cluster were found 
to decrease as a power law with time, with the velocity 
fluctuations being much smaller than the centre of mass 
velocities. However, when these fluctuations became of 
order 5, then intra cluster collisions became mostly elas- 
tic and the clusters start to fragment. This emergent 
phenomena was quantified by defining an effective frag- 
mentation rate for a cluster. The fragmentation rate 
was seen to rise sharply at some cluster size dependent 
time with the crossover time increasing with decreasing 5. 
Once fragmentation sets in, the cluster size distribution 
{N{m,t)) changes drastically from that of the sticky gas 
(ro — 0, S — 0). However, the scaling of {N{m,t)) could 
be related to that of the velocity fluctuations a{m, t) 
which was completely characterised by two independent 
exponents xi and X2- It was also observed that the total 
energy of the system as well as clusters continue to decay 
as t~'^^^, showing no signature of the structural changes 
in the clusters. 

We believe that many of these results (qualitative) will 
be carried over to higher dimensions. In two dimensions, 
for coarse grained velocities, it was shown ^] that the 
velocity fluctuations scale differently from the typical ve- 
locity. Hence, we expect a crossover when these fluctu- 
ations become comparable to S, and thus fragmentation 
to be relevant for two dimensions too, and consequently 
a regime where coarsening is fluctuation dominated. It 
would be interesting to verify it numerically. 

The velocity scale 6 is relevant and not just a computa- 
tional tool. Experimentally, r(v) approaches 1 when the 
relative velocity v tends to zero, i.e., r{v) ~ 1 — (v/S)^ + 
. . ., v/S 1. The exponent x takes a variety of values. 
Within the framework of viscoelastic theory, x — 1/5- 
Systems with x < 1 cannot be studied using event driven 
molecular dynamics simulations as inelastic collapse pre- 
vents the simulation from proceeding forward. It would 
be interesting to use conventional molecular dynamics 
simulations to verify whether the observations of this pa- 
per is valid for % < 1. 

Another question of interest is the construction of lat- 
tice models that reproduce the coarse grained behaviour 
of the granular gas. Such models are not only computa- 
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tionally much faster, but also may be the first steps to- 
wards building effective field theories for the system. In a 
recent paper a stochastic lattice model was studied 
which reproduced all features of the sticky gas. It would 
be interesting to see whether fragmentation can be in- 
corporated into this lattice model such that the coarse 
grained behaviour seen in this paper is reproduced. 



Academy of Sciences under the Grant No. 2011001. A 
part of this work (MS) was carried out at Centre for 
Nonlinear Studies, Hong Kong Baptist University. 



Acknowledgments 



MS was partially supported by the International Young 
Scientist Fellowship of Institute of Physics, Chinese 



[1] G. F. Carnevale, Y. Pomeau, and W. R. Young, Phys. 

Rev. Lett. 64, 2913 (1990). 
[2] S. Mcnamara and W. R. Young, Phys. Fluids. A 4, 496 

(1992) . 

[3] I. Goldhirsch and G. Zanetti, Phys. Rev. Lett. 70, 1619 

(1993) . 

[4] L. Frachebourg, Phys. Rev. Lett. 82, 1502 (1999). 

[5] L. Frachebourg, P. A. Martin, and J. Piasecki, Physica 

A 279, 69 (2000). 
[6] E. Ben-Nairn, S. Chen, G. Doolen, and S. Redner, Phys. 

Rev. Lett. 83, 4069 (1999). 
[7] X. Nie, E. Ben-Nairn, and S. Chen, Phys. Rev. Lett. 89, 

204301 (2002). 

[8] E. Trizac and A. Barrat, Eur. Phys. J E 3, 291 (2000). 
[9] E. Trizac and P. L. Krapivsky, Phys. Rev. Lett. 91, 
218302 (2003). 
[10] S. Das and S. Puri, Phys. Rev. E 38, 011302 (2003). 
[11] M. Shinde, D. Das, and R. Rajesh, Phys. Rev. Lett. 99, 
234505 (2007). 

[12] M. Shinde, D. Das, and R. Rajesh, Phys. Rev. E. 79, 

021303 (2009). 
[13] P. Haff, J. Fluid Mech. 134, 401 (1983). 
[14] C. C. Maa,g, N. Isert, G. Maret, and C. Aegeter, Phys. 

Rev. Lett. 74, 248001 (2008). 



[15] S. Kida, J. Fluid Mech. 93, 337 (1979). 

[16] G. Porod, "Small-angle x-ray scattering," (Academic, 

London, 1982). 
[17] C. V. Raman, Phys. Rev. 12, 442 (1918). 
[18] L. Labous, A. Rosato, and R. Dave, Phys. Rev. E 56, 

5717 (1997). 

[19] E. Falcon, C. Laroche, S. Fauve, and C. Coste, Euro. 
Phys. J. B 3, 45 (1998). 

[20] N. V. Brilliantov and T. Poschel, Kinetic Theory of Gran- 
ular Gases (Oxford University Press, 2004). 

[21] C. Huepe and M. Aldana, Phys. Rev. Lett. 92, 168701 
(2004). 

[22] N. G. V. Kampen, Stochastic Processes in Physics and 

Ghemistry (Elsevier, Netherlands, 2007). 
[23] I. Pagonabarraga and E. Trizac, Lect. Notes Phys. 624, 

163 (2003). 

[24] R. C. Hidalgo and I. Pagonabarraga, Europhys. Lett. 77, 
64001 (2007). 

[25] D. C. Rapaport, The art of molecular dynamics simula- 
tions (Cambridge University Press, Cambridge, 2004). 

[26] S. Dey, D. Das, and R. Rajesh, Euro. Phys. Lett. 93, 
44001 (2011). 



